Dynamical collapse in a degenerate binary fermion 
mixture using a hydrodynamic model 

Sadhan K. Adhikari| 

Institute de Fi'sica Teorica, UNESP - Sao Paulo State University, 01.405-900 Sao 
Paulo, Sao Paulo, Brazil 

Abstract. 

We use a time-dependent dynamical hydrodynamic model to study a collapse in 
a degenerate fermion-fermion mixture (DFFM) of different atoms. Due to a strong 
Pauli-blocking repulsion among identical spin-polarized fermions at short distances 
there cannot be a collapse for repulsive interspecies fermion-fermion interaction. 
However, there can be a collapse for a sufficiently attractive interspecies fermion- 
fermion interaction in a DFFM of different atoms. Using a variational analysis and 
numerical solution of the hydrodynamic model we study different aspects of collapse 
in such a DFFM initiated by a jump in the interspecies fermion-fermion interaction 
(scattering length) to a large negative (attractive) value using a Feshbach resonance. 
Suggestion for experiments of collapse in a DFFM of distinct atoms is made. 
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1. Introduction 

Recent successful observation of degenerate boson-fermion mixture (DBFM) and 
fermion-fermion mixture (DFFM) of trapped alkali-metal atoms by different 
experimental groups [1-7] has initiated the intensive experimental studies of different 
novel phenomena [8-10]. It has been possible to achieve a degenerate Fermi gas (DFG) 
by sympathetic cooling in the presence of a second boson or fermion component, as 
there cannot be an effective evaporative coohng [1] of a single-component DFG due to 
a strong Pauli-blocking repulsion at low temperature among spin-polarized fermions. 
Among these experiments on a DFG, apart from the study of a DBFM in '^'^Li [6], 
23Na-^Li [7] and ^^Rb-^^K [9-12], there have been studies of a DFFM in 40k-40k* [1] 
and ^Li-^Li* [2-5] systems, where * denotes a distinct hypcrfine state. More recently the 
formation of a Bardeen-Cooper-Schreiffer (BCS) condensate of fermionic ®Li atoms in 
^^Na-^Li [13] and ^Li-^Li* [4] mixtures has been observed experimentally. The collapse 
in a DBFM of ^^Rb-^°K atoms has been observed and studied by Modugno et al. [9] and 
more recently by Ospelkaus et al. [11]. Recently, experiments on controlled collapse on 
^^Rb-^°K have been accomplished [14] . 

A collapse in a Bose-Einstein condensate (BEG) takes place due to an attractive 
atomic interaction [15, 16]. A study of controlled collapse and explosion has been 
performed by Donley et al. [16] on an attractive *^^Rb BEG, where they manipulated 
the inter-atomic interaction by varying a background magnetic field exploiting a nearby 
Feshbach resonance [17]. There have been many theoretical [18, 19] studies to describe 
different features of this experiment [16]. More recently, there have been experimental 
studies on collapse in a DBFM of ^'''Rb-^^K by two different groups [9, 11, 14] as well 
as related theoretical investigations [20,21]. As the interaction in a pure DFG at short 
distances is repulsive due to Pauh-blocking, there cannot be a collapse in it. The 
Pauli repulsion is responsible for the stability of a neutron star against a (gravitational) 
collapse. A collapse is possible in a DBFM in the presence of a sufficiently strong boson- 
fermion attraction which can overcome the Pauli repulsion among identical fermions 
[9,11,20]. 

In this paper we study the collapse in a DFFM for a sufficiently attractive 
interspecies fermion-fermion interaction which can overcome the Pauli repulsion. 
However, there is already experimental evidence and theoretical conjecture that a Fermi 
gas in two (spin) hyperfine states of the same atom is much more stable [22] than 
expected on the basis of a scattering length approach [5, 23] and there is no collapse in 
such a system. A similar conclusion follows from an examination of the compressibility 
of such a system [24] . A strongly attractive DFFM in two (spin) hyperfine states exhibits 
universal behavior and should be mechanically stable as a consequence of the quantum- 
mechanical requirement of unitarity. This requirement limits the maximum attractive 
force for such a DFFM to a value smaller than that of the outward Fermi pressure 
due to Pauh repulsion [22] . It has been demonstrated that a two-component DFFM in 
different (spin) hyperfine states is stable against collapse [25] , whereas a multicomponent 
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degenerate fermion mixture [22, 25] or a DFFM of different atoms [26] could undergo 
collapse. Hence, by taking a DFFM of two different fermionic atoms of different atomic 
mass one can avoid the problem [22, 25] of a possible suppression of collapse. Thus 
one could study collapse in a DFFM in the same manner as in a DBFM. The second 
component of fermions will then only aid in inducing an attraction among the fermions of 
the first component responsible for collapse without suppressing the collapse. Although, 
the past experiment [22] on two-component cold Fermi gas used two (spin) hyperfine 
states of the same atom, experiments can be realized with distinct atoms and one can 
look for collapse in such a system. One such system is the ^Li-^°K mixture: both ^Li [4] 
and [1] have been trapped and studied in laboratory. 

Here we use a coupled time-dependent mean-field hydrodynamic model which is 
inspired by the success of a similar model used by the present author in the investigation 
of a fermionic collapse [20] and bright [27] and dark [28] solitons in a DBFM as well 
as of mixing-demixing [29] and black solitons [30] in a DFFM. The conclusions of the 
study on bright soliton [27] are in agreement with a microscopic study [31] and the noted 
survival of collapse in the numerical study [20] has been experimentally substantiated 
later in a DBFM of ^°K-^^Rb [11]. A very similar model has been used by Jezek et 
al. [32] in a successful description of vortex states in a DBFM. Although, a mean-field 
model is simple to use and leads to a proper prediction of probability density of the 
fermionic system, many true quantum effects are lost in this simplified model, e. g., it 
cannot predict the suppression of collapse of a DFFM in two different (spin) hyperfine 
states [22,25] as discussed in the last paragraph. We recall many true quantum effects 
are also lost [33] in the mean-field Gross-Pitaevskii equation for trapped bosons. 

In our study on collapse in a DFFM we shall consider a strong attraction among 
fermions which will naturally lead to molecule formation and not to a BCS state. A BCS 
state is usually formed for a weak attraction among identical fermions. The possibility 
of molecule formation by three-body recombination is explicitly included in our model 
by an absorptive nonlinear term. Apart from the direct experimental interest in trapped 
cold atoms, a study of strongly interacting Fermi gases and their possible collapse is also 
relevant [22] in condensed matter physics (superconductivity), nuclear physics (nuclear 
matter), high energy physics (effective theories of strong interaction), and astrophysics 
(compact stellar objects), which makes the present study of greater interest. 

The collapse in a DFFM is first studied by a variational analysis of the present 
model, which is later substantiated by a complete numerical solution using the Crank- 
Nicholson scheme [34]. During a collapse and an explosion of the DFFM, the system 
loses atoms as in a collapsing and exploding BEC [16]. The loss of atoms is accounted for 
by three-body recombination involving two types of fermions. We study the sensitivity 
of our results on the three-body recombination loss rates. We also study the quasi- 
periodic oscillation of the sizes of the DFFM undergoing a collapse. The collapse is to 
be initiated by jumping the interspecies scattering length to a large negative (attractive) 
value near a fermion-fermion Feshbach resonance [35]. However, the collapse starts after 
a time delay upon this jump and we study the variation of this time delay with the final 



Collapse in a fermion-fermion mixture 



4 



scattering length. This variation has a behavior similar to that observed in the bosonic 
case [16]. The collapsing DFFM is found to execute a quasi-periodic oscillation with 
a frequency approximately equal to twice the frequency of the harmonic trap as in a 
BEC [16]. 

Previously, in addition to the study of a collapse in a pure BEC [18, 19], we also 
investigated [36] the collapse in a two-species BEC initiated by an interspecies attraction. 
The predicted collapse in a two-species BEC for intra-species repulsion and interspecies 
attraction [36] is similar to that in a DBFM of ^^Rb-^^^K studied before [9,11,20] and 
that in a DFFM studied in this paper. 

In section 2 we present the coupled hydrodynamic model for a DFFM which we 
apply to predict and study a collapse. In this section we also present a variational 
analysis based on this model which substantiates the collapse in a DFFM for a 
sufficiently attractive interspecies fermion-fermion attraction. In section 3 we present 
results of numerical simulation of our study on collapse. Finally, in section 4 we present 
a brief summary of our investigation. 

2. Coupled Hydrodynamic Model for a Fermion-Fermion Mixture 

2.1. Model 

A mean-field-hydro dynamic Lagrangian for a DFG has been used successfully in the 
study of a DBFM [27,28,32] which we shall use in the present investigation. The virtue 
of the hydrodynamic model for a DFG over a microscopic description is its simplicity 
and good predictive power. To develop a set of practical time-dependent hydrodynamic 
equations for a DFFM, we consider the following Lagrangian density [27, 28] 

2 ^ dt ' dt J 

i=l,2 ^ / 

j=l \ J / 
+ 9i2nin2 - ih (^31^1^^ + ^32^^712) , (2.1) 

where rrij is the mass of component j(= 1,2), Aj = ^^(67r^)^/^/(2mj), \I'j is a complex 
probability amplitude, nj = is a real probability density and Nj = J drnj{r) 

the number. Here the interspecies coupling is gu = 2TTh'^ai2/mR with the reduced 
mass rrifi = 71111712/ {mi + 777.2), and ai2 is the interspecies scattering length. The 
spherically-symmetric potential is taken as V,(r) = |(3mj)co'^r^ where uj is the radial 
(r) frequency. The interaction between intra-species fermions in spin-polarized state 
is highly suppressed due to Pauli blocking and has been neglected in fl2.H) and will 
be neglected throughout. The kinetic energy terms in this equation are derived from 
a hydrodynamic equation for the fermions [37]. The kinetic energy terms contribute 
little to this problem compared to the dominating Pauli blocking term 3y4j|nj|^/^/5 
in 1)2.11) . However, its inclusion leads to a smooth solution for the probability density 
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everywhere [27]. The Lagrangian density of each fermion component in ()2.H) is identical 
to that used in Refs. [27,28]. The last two terms in 1)2.11) correspond to three-body 
recombination due to the following reactions, respectively F1+F2+F2 — (FiF2)+F2, and 
F1 + F2 + F1 {F1F2) +Fi, where (F1F2) is a composite structure (resonance/molecule) 
of fermions Fi and F2 and K^i and K32 are the corresponding three-body loss rates. The 
contribution to the Lagrangian density of the recombination reactions is proportional 
to the density of the participating fermions. Here we neglected two-body loss. 

The dynamical equations for the system are just the usual Euler-Lagrange (EL) 
equations with the Lagrangian density ()2.1|) [38] 

d dC dC 



d dC 



dtd 



^ dXk 

dt k=l " ^ dxf, 

where Xk,k = 1, 2, 3 are the three space components, and j = 1,2 refer to the fermion 
components. Consequently, the following EL equations of motion are derived: 



6mi 

ifi (Kaing + 2i^32nin2) 



^1 = 



2/3 



0. 



+ 9X2^2 



(2.3) 



ih^ - 4^ + V2{j) + A2\n2fl^ + gx2nx 



dt 6m2 

- ih (2Ksinin2 + 



*2 = 0. 



(2.4) 



In the spherically-symmetric state the fermion density has the form \E'j(r;t) = 
tpj{r]t). Now transforming to dimensionless variables defined hy x = ^/2r/l, r = tco, 
I = \/ h/ [muj), m = 3mi = 3m2 and 



Pj[x;T) 



47rP 



we obtain from (lOl and 

d 



ipj{r;t), 



(2.5) 



d 


2 


dx^ 




2 


02 

I 




X 






4 


02 




X 





X 

T 



4/3 



01 (x; r) = 0, 





2 




01 




02 


X 




X 



(2.6) 



d 92 



— I 



dr dx'^ 



X 

T 



22 



X 



4/3 



+ QV2M21 


01 


2 

- i2UMl2M2l 


01 


2 


02 




X 


X 




X 



ii^2M, 





4 


01 




X 





52 la;; rj 



(2.7) 



Collapse in a fermion-fermion mixture 



6 



where ^^Jj = 3(37riV,/2)2/3, ^f^^ = N2ai2/l, M21 = N^ai^/l, 62 = K^2/ {'iT^^al^l^uj), 
and = ir3i/(27r^a^2^^^)- ^^e non-absorptive case without any loss of atoms due 
to three-body recombination ^31 = ^32 = 0, the normahzation of the wave-function 
components is given by r)p = l,j = 1,2. In the absorptive case ^31 7^ 

and (^32 7^ it is possible to have loss of atoms due to three-body recombination and 
the normalization reduces with time due to loss of atoms. 

We solve the coupled hydrodynamic equations ()2.6p and ()2.7p numerically using a 
time-iteration method based on the Crank-Nicholson discretization scheme elaborated 
in Refs. [34,39]. We discretize the hydrodynamic equations using time step 0.00025 
and space step 0.025 spanning x from to 25. This domain of space was sufficient 
to encompass the entire fermion function during a collapse and explosion and obtain 
convergent solution for the total number of fermions. First we solve ()2.6|1 and ()2.7|1 
with ^31 = ^32 = to find an initial stationary state of the DFFM. It is true that the 
three-body loss, taken care of by terms ^31 and ^32, is always present in the system, 
its effect is small leading to at best a small loss rate in atoms except for very large 
negative values of 012- Hence, for the consideration of the initial state, we could as well 
neglect three-body loss. (This is why a Gross-Pitaevskii mean-field equation without 
three-body loss has been successfully used to study many features of a repulsive trapped 
BEC [33].) However, in the final state, when ai2 is suddenly turned to a large negative 
value by a Feshbach resonance, the three-body loss has a dramatic effect responsible for 
a proper description of a collapse and explosion with a very rapid loss of atoms in a 
short interval of time (see figure 3). 

In our numerical investigation we take / = 1 /im and consider the equal-mass 
fermions with the mass of ^"K corresponding to a radial frequency «i 27r x 83 Hz. The 
present simplified mean-field model cannot predict the suppression of collapse [22,25] of a 
DFFM in two hyperfine states which is a true quantum many-body effect. Nevertheless, 
it leads to a proper description of collapse dynamics of a DFFM of two distinct atoms. 
The use of equal-mass fermions only keeps the algebra simple specially in section 2.2, but 
otherwise has no effect on the general qualitative dynamics studied in this paper. In this 
study, the unit of time is l/cj ^ 2 ms, and unit of length 1/ \/2 ^ 0.7 yum. Actually, any 
two different fermionic atoms can be used in experiment, a proper quantitative treatment 
of which will require the use of different mass factors in the dynamical equations. 



2.2. Variational Analysis 



To understand how the stationary states of a DFFM are formed, we employ a variational 
method for the solution of ()2.6|) and ()2.7p in the symmetric case Ni = N2 = N, while 

(pi/x = (f)2/x = ip satisfies 



d 



2 d 



X 



1 4/3 



+ 9 



0, 



(2.^ 



dr dx"^ 

where g = Q\/2Nai2/l and /i = 3(37rA^/2)^/'^ [40]. Here we have set the absorptive terms 
to zero for stationary states. We consider the following trial Gaussian wave function for 
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ip{x, t) = A{t) exp 



;2.9) 



where A{t), R{t), (3{t), and a{t) are the normahzation, width, chirp, and phase, 
respectively. The normalization condition dxx^ip'^{x,t) = 1 sets A(t) = 
[7r^/^i?^/^(i(:)/2]^^. The Lagrangian density for generating ()2.8|) is [40] 



dip 



dx 



X 



(2.10) 



where the overhead dot represents time derivative. The trial wave function ()2.9|) is 
substituted in the Lagrangian density and the effective Lagrangian L^g is calculated 
via Lgg = J C{{p)(Px : 



^eff 



7i^/^A\t)R'{t) 



9V3 A^/\t) 2a{t) 
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The generalized Lagrange equations for this effective Lagrangian given by [38 



hit) - 

2^^ ' 2V2R\t) 
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with 7(t) representing a{t), A{t),i3{t), and R{t) are written explicitly as 
^3/2^2_^3 ^ constant = An, 



R = 2Rp, 



i?2 5^5 R' 



gA^ 18v^/iA4/3 
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(2.12) 
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(2.16) 



where the time dependence of different observable is suppressed. Eliminating a between 
flTTI|l and ^TT^i one obtains 



2$ = —,-4f3' 



gA' ^ uVSfiA'/' ^ 



(2.17) 



V2R^ 25^5 R^ 

From ()2.15j) and ()2.17j) we get the following second-order differential equation for 
the evolution of the width R 



d^R 
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Ag 12^42/3^3 



dt^ i?3 257ri/3v/5i?3 



R. 



(2.18) 
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2 4^ 
/?2 + 3 



1 6/i42/3v/3 R 



+ 



(2.19) 



257ri/375i?2 

The quantity in the square brackets of ()2.19|) is the effective potential U{R) of the 
equation of motion: 



6/i4^/3v/3 ^ 
^ 257ri/3^/?2 + 2 ■ 



(2.20) 



Small oscillation of a stationary state around a stable configuration is possible when 
there is a minimum in this effective potential determined by a zero of ()2.18|) . This 
condition yields the variational width from which the variational solution for the wave 
function is obtained via ()2.9p . 

In figure 1 we plot the effective potential U{R) of ()2.2()j) for different ai2, N = 1000 
and / = 1 jum. For positive (repulsive) = 100 nm, U{R) leads to a confining well with 
a minimum at R = Rq = 5.3, so that one could have a stable DFG of width Rq = 5.3. 
The variational profile for this function is 



3/2 



exp 



X 



2 1 



2RI 



f2.211 



In figure 1, as turns gradually negative (attractive), the infinite wall near i? = of 
U{R) is gradually lowered and for a sufficiently attractive scattering length ai2 ~ —100 
nm, this wall disappears completely and one has the possibility of collapse for ai2 < —100 
nm. The minimum in U{R) first becomes a point of inflection for ai2 ~ —100 nm and 
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Q [ ' .Stable bosons 
-0.1 -0.05 

ai2/l 

Figure 2. The phase diagram for collapse where we plot logio(A^) vs. 012// for 
fermion-fermion mixture (full line) and bosons (dotted line). The line separates the 
regions where collapse is present and absent. 

then disappears. The profile of the effective potential in figure 1 is similar to the same 
in the bosonic case [33,40]. 

Next we show the phase diagram for collapse for Ni = N2 = N using the variational 
approach in figure 2, where we plot logio(A^) vs. ai2/l. The line separates the plot in 
two regions. In the upper half of the plot collapse is possible and in the lower half 
we have stable configurations. The phase diagram of figure 2 is quite similar to that 
obtained in Ref. [23] in a study of the stability of a DFFM. In case of bosons in a 
spherically-symmetric trap the line of stability is given by Na/l = —0.575 [33] and is 
also shown in figure 2. As expected, for a fixed |a//| a much larger number of fermions 
can be accommodated in a stable state. 

3. Numerical Results 

In this section we present results on collapse from a numerical solution of ()2.6p and 
(jTTj). After some experimentation we take in the initial DFFM A^i = 1000, N2 = 2000, 
and ai2 = 100 nm, so that 012// = 0.1. This corresponds to nonlinearities A/11 = 843, 
A/22 = 1338, A/12 = 200 and A/'2i = 100. The collapse dynamics is sensitive to the loss 



Collapse in a fermion-fermion mixture 



10 



1000 



500 







2 
1 



K = 1 0'^^cm^s"^ 




10 cm S 



1 0'^^cm^s"^ 







10 20 30 
t(ms) 



40 



2000 



CM 



1000 







50 



Figure 3. The evolution of fermion numbers Nj{t) of the two components vs. time 
during a coUapse initiated by a jump in scattering length ai2 from 100 nm to —200 nm 
in a DFFM of iVi 1000 and N2 = 2000 for three-body loss rates K = IQ-^e cm*^s-\ 
10^^^ cm^s^^, and 10^^'' cm^s^^. The dotted (blue) curves refer to fermion 2 and the 
solid (red) curves to fermion 1. The curves are labeled by their respective K values. 



rates K31 and -ft'32. As these loss rates are not experimentally known, in the present 
simulation we take them to be equal: K = K31 = K32, and consider different values of 
K. 

Now we consider the collapse of fermions initiated by a sudden jump in the fermion- 
fermion scattering length from 012 = 100 nm to —200 nm which can be implemented 
near a fermion-fermion Feshbach resonance, observed in fermionic systems [35]. This 
resonance should enable an experimental control of the interspecies interaction [17] 
and hence can be used to increase the attractive force between interspecies fermions by 
varying a background magnetic field, which in turn increases the attractive nonlinearities 
6-\/2A/i2 and 6-\/2A/2i in ()2.fi|) and ()2.7j) . If these attractive nonlinear terms can overcome 
the repulsive nonlinearities in these equations it is possible to have a collapse of fermions. 

Due to the three-body loss terms in ()2.(j|l and ()2.7|1 the number of fermions decay 
with time. When the net nonlinear attraction in these equations is small there is 
a smooth and steady decay of number of atoms. However, when the net nonlinear 
attraction is jumped to a large value, the steady decay of number of atoms develops 
into a violent decay called collapse. When this happens, the DFFM loses a significant 
fraction of atoms in a small interval of time (milliseconds) after which a remnant DFG 
with a reasonably constant number of atoms is formed. Also, during and immediately 
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Figure 4. The fermion density nj{r) a,t t — 0,20 ms before and during the collapse 
exhibited in figure 3 for K — 10~^^ cm^ s~^. The density calculated from the 
variational profile of the wave function H2.21|l for Rq — 5.3 and N = 1000 is also 
shown, 



after collapse, the fermion density function becomes unsmooth and spiky in nature 
as opposed to a reasonably smooth function in the case of a steady decay. This also 
happened in the collapse of a BEC [16]. 

We study the evolution of fermion numbers in the DFFM from time t = to t = 50 
ms after a sudden jump in the scattering length from ai2 = 100 nm to —200 nm at t = 0. 
In agreement with the variational analysis of last section we find that this jump in the 
scattering length leads to collapse. The evolution of fermion numbers after the jump in 
scattering length au depends on the value of the three-body loss rate K. We study the 
sensitivity of the result on K by performing the calculation for different loss rates. In 
figure 3 we plot the evolution of the fermion numbers for loss rates: K = 10^^^ cm^s~^, 
10^^^ cm^s~^, and 10~^^ cm^s~^. With the increase of K, the decay rate increases, 
although the results for different K are qualitatively similar. We see in figure 3 that, in 
all cases, the number of fermions decays rapidly and attain an approximately constant 
(remnant) value in less than 50 ms as in the case of bosons [16, 18]. We used a space 
step of 0.025 in the numerical solution of ()2.6|1 and ()2.7|1 and found that this step size 
was sufficient to reach a converged result even during collapse. In our previous study 
on the collapse of a BEC of ^^Rb [18, 20] we found that even a much larger step size 
of 0.1 led to converged result for the number of atoms in the remnant in quantitative 
agreement with experiment [16]. This gives assurance on the reliability of the present 
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Figure 5. Evolution of central density nj{0) of fermion j = (a) 1 and (b) 2 during 
the collapse exhibited in figure 3 for = 10~^^ cm^s~^. 



calculation. 

In figure 4 we plot the fermion probability densities at times t = and t = 20 ms. 
A close look at figure 4 reveals that before collapse at t = the fermion densities are 
smooth. We have also plotted here the density corresponding to the variational profile 
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()2.21|1 for Ro = 5.3 for fermion number N = 1000. Although the ranges of the exact 
and variational densities agree with each other, the functions do not agree well. This is 
understandable as the exact profile of the DFG is very different from the Gaussian trial 
function used in variational calculation. Although the fermion densities at t = are 
smooth, the fermion densities during and after collapse have an entirely different profile. 
As expected the densities are highly peaked in the central (r = 0) region and develop 
spikes. Near r = the densities could be two to three orders of magnitude larger than 
those for larger r values (see figure 5). However, they extend over a large distance too. 
The final spiky function indicates the collapse, in contrast to a smooth final function 
corresponding to a steady loss of atoms. The collapse is a quick process lasting at most 
a few tens of milliseconds when a significant fraction of atoms are lost. For example, in 
figure 3 for K = 10"^"^ cm^s~^, the collapse lasts for the first 25 ms when most of the 
atoms are lost. After this interval the rate of loss of atoms is reduced and remnant a 
DFG with a roughly constant number of atoms are formed. 

To confirm further the collapse in figures 3 and 4 for i^' = 10~^^ cm^s"^, we plot 
in figure 5 the evolution of the central probability density nj{0) of fermion j during 
collapse. We note a very strong fluctuation of a very large central density reminiscent 
of collapse in both components. The central density is three orders of magnitude larger 
than the equilibrium density in figure 4. Similar fiuctuations were noted in the collapse 
of a pure BEG [19] as well as a DBFM [20]. Such a strong fluctuation of the central 
density could not be due to a weak evaporation of the DFFM due to recombination. 

From flgure 3 we flnd that the number of fermions remains practically constant 
during the flrst 4 ms or so after jumping the scattering length from 100 nm to —200 
nm indicating that the collapse starts only after this interval of time. This is conflrmed 
from the plot of central densities in flgure 5 where we see that very large values of 
central density also appear after an interval of time called "time to collapse" . A similar 
phenomenon was also observed in the collapse of bosons [16]. Next we study an evolution 
of this time to collapse ( ''^collapse ) ^i^^ changing initial (^initial) iinaA (o^Qjjg^pgg) 
scattering lengths. This is shown in flgure 6 for two values of ^initial' '^here we plot 
^collapse "collapse = ^^00, N2 = 2000 and K = lO^^e ^.^e^-i^ rj.^^ ^-^^^ 

collapse is large for a small jump in the scattering length and reduces when the jump 
in the scattering length is increased, as also observed in the case of bosons [16]. 

One interesting aspect of flgure 3 is the appearance of a revival of collapse. The 
fermion number after the primary collapse remains approximately constant for an 
interval of time and then again reduces abruptly. This revival of collapse takes place 
several times. A similar revival of collapse was noted in the fermion component in a 
numerical simulation in a DBFM [20] and was conflrmed later in an experiment [11] on 
the ^''Rb-^'^K DBFM. To study the revival of collapse further in a DFFM we considered 
a different jump in the scattering length. For the same initial state of flgure 3 we now 
consider a jump in ai2 from 100 nm to —300 nm and the dynamics is reported in flgure 
7 for different K values. We flnd that the revival of collapse has practically disappeared 
in this case. If the collapse is initiated by a small jump in the scattering length, the 
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Figure 6. The evolution of time to collapse ^collapse ^^^^^^ scattering length 
"coUapse ^initial = ^ and 100 nm for a DFFM with A^i 1000, iVa 2000 and 
K — 10~^^ cm^s~^. The curves are labeled by their respective ajnitial "^^-lues. 

initial collapse is less violent. However, after this initial milder collapse the DFFM 
cannot reach an equilibrium state and it remains large and attractive. Consequently, 
the DFFM undergoes further collapse(s). On the other hand, if the collapse is initiated 
by a large jump in the scattering length, the initial collapse is very violent through 
which the DFFM gets rid of a very large number of atoms. Consequently, the DFFM 
reaches reasonably small and cold remnant states which do not further undergo collapse 
and one has one primary collapse. This is clear from figures 3 and 7. In figure 3 after 
the first collapse the DFFM loses a smaller percentage of atoms whereas in figure 7 a 
large percentage of atoms are lost after the primary collapse. 

It was found in the experiment on collapse on a BEC [16] that during collapse 
the root mean square (rms) sizes of the condensate execute periodic oscillation with 
approximately twice the frequency of the trap. A breathing oscillation of same frequency 
was found in a BEC when a small perturbation was applied [34]. In dimensionless unit, 
the angular frequency of the trap is = 1, corresponding to a (linear) frequency of 
l/(27r). The observed frequency of oscillation of rms size was l/vr [16] — twice the trap 
frequency. In actual time unit the frequency of oscillation of rms sizes during collapse 
corresponds to 1/ (27r) ms^^ ^ 0.16 ms~^. We also investigated if such oscillation existed 
in the present case in the DFFM during the collapse. In figure 8 we plot the rms radii 
of the two components of the collapsing DFFM and find that they also execute quasi- 
periodic oscillation. The calculated frequency from figure 8 is 0.145 ms~^ close to that 
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Figure 7. The evolution of fermion numbers Nj{t) during collapse initiated by a 
jump in scattering length ai2 from 100 nm to —300 nm for K = 10~^^ cm^s~^, 10~^^ 
cm^s~^, and 10~^^ cm^s"^. The dotted (blue) curves refer to fermion 2 and the solid 
(red) curves to fermion 1. The curves are labeled by their respective K values. 

found in the case of bosons, e.g. 0.16 ms~^. The difference could be due to the coupled 
nature of the hydrodynamic equations as well as the very large nonlinearity for fermions. 

4. Summary 

We suggested a coupled set of time- dependent hydrodynamic equations for a trapped 
DFFM including the effect of three-body recombination. The present time- dependent 
formulation permits us to study the non-equilibrium dynamics of a DFFM. Using a 
variational analysis as well as a numerical solution of our model, we study, for an 
attractive inter-species fermion-fermion interaction, the collapse in a DFFM composed 
of two types of nonidentical atoms. The collapse of a DFFM of two different atoms 
can be realized experimentally by jumping the inter-species scattering length to a large 
negative value by exploiting a fermion-fermion Feshbach resonance [35]. The collapse 
dynamics is strongly dependent on the three-body loss rate K and we present results 
for different loss rates. We note the possibility of a revival of collapse in a DFFM as in a 
previous simulation [20] on a DBFM, confirmed later in an experiment [11] on ^^Rb-'^^K. 
We find that a revival of collapse in a DFFM takes place for a moderate jump in the 
interspecies scattering length which disappears for a larger jump. We also study the 
quasi-periodic oscillation of the DFFM with approximately twice the trap frequency 
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Figure 8. The rms radii of the two components vs. time before and during cohapse 
exhibited in figure 3 for A^i = 1000 and N2 = 2000 initiated by a jump in scattering 
length ai2 from 100 nm to —200 nm for K = 10^^^ cm^s^^. The dotted (blue) curves 
refer to fermion 2 and the solid (red) curves to fermion 1. 

during collapse and explosion. 

A proper treatment of a DFFM should be performed using a fully antisymmetrized 
many-body Slater determinant wave function [41] as in the case of atomic and molecular 
scattering involving many electrons [42]. However, in view of the success of the 
hydrodynamic model in a description of a collapse [20], the formation of bright [27] 
and dark [28] solitons, and vortex states [32] in a DBFM, we do not believe that the 
present study on the collapse in a DFFM to be so peculiar as to have no general validity. 
The present study on collapse in a DFFM of nonidentical atoms can be verified in future 
experiments, which can really validate the present hydrodynamic model and the related 
numerical study. 
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